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ABSTRACT 

We investigate the origin of the soft X-ray excess emission from narrow-line Seyfert 1 galaxies 
Akn 564 and Mrk 1044 using XMM-Newton observations. We find clear evidence for time delays 
between the soft and hard X-ray emission from Akn 564 based on a ~ 100 ks long observation. 
The variations in the 4-10 keV band lag behind that in the 0.2 - 0.5 keV band by 1768 ± 122 s. 
The full band power density spectrum (PDS) of Akn 564 has a break at ~ 1.2 x 10~ 3 Hz with 
power-law indices of ~ 1 and ~ 3 below and above the break. The hard (3 — 10 keV) band PDS 
is stronger and flatter than that in the soft (0.2 — 0.5 keV) band. Based on a short observation 
of Mrk 1044, we find no correlation between the 0.2 — 0.3 keV and 5 — 10 keV bands at zero lag. 
These observations imply that the soft excess is not the reprocessed hard X-ray emission. The 
high resolution spectrum of Akn 564 obtained with the reflection grating spectrometer (RGS) 
shows evidence for a highly ionized and another weakly ionized warm absorber medium. The 
smeared wind and blurred ionized reflection models do not describe the EPIC-pn data adequately. 
The spectrum is consistent with a complex model consisting of optically thick Comptonization in 
a cool plasma for the soft excess and a steep power-law, modified by two warm absorber media 
as inferred from the RGS data and the foreground Galactic absorption. The smeared wind and 
optically thick Comptonization models both describe the spectrum of Mrk 1044 satisfactorily, 
but the ionized reflection model requires extreme parameters. The data suggest two component 
corona - a cool, optically thick corona for the soft excess and a hot corona for the power-law 
component. The existence of a break in the soft band PDS suggests a compact cool corona that 
can either be an ionized surface of the inner disk or an inner optically thick region coupled to a 
truncated disk. The steep power-law component is likely arising from an extended region. 

Subject headings: accretion, accretion disks - galaxies: active - X-rays: galaxies 



1. Introduction 

A significant fraction of type 1 active galactic 
nuclei (AGN) show 'soft X-ray excess emission' 
above a power-law continuum, usually identified 
as the steepening of the X-ray continuum below 
~ 2 keV. This soft excess emission was first ob- 
served by HEAO-1 (Pravdo et al. 1981), and EX- 
OSAT (Arnaud et al. 1985; Singh et al. 1985). 
Boiler, Brandt, & Fink (1996) showed that AGN 
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with the steepest soft X-ray spectra tend to lie 
at the lower end of the H/3 line width distribu- 
tion. These AGN with FWHM H(3 £ 2000 km s" 1 
are classified as the narrow-line Seyfert 1 galaxies 
(NLS1; Osterbrock & Pogge 1985), and are dis- 
tinguished from the bulk of the Seyfert 1 galaxies. 
In addition to the strong soft X-ray excess below 
2 keV and narrower permitted lines, NLS1 also 
show a number of extreme properties e.g., steep 
2 — 10 keV power-law continuum, extreme X-ray 
variability and strong Fe II emission. These prop- 
erties are all related and one suggestion is that 
NLS1 are the supermassive black hole analogue of 
stellar mass black hole X-ray binaries (BHB) in 
their 'high/soft' state. 
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However, the origin of the soft X-ray excess 
emission has remained a major problem in AGN 
research over the past two decades. The soft ex- 
cess emission is a smooth continuum component 
rather than a blend of emission/ absorption fea- 
tures as revealed by Chandra and XMM-Newton 
grating observations (e.g., Turner et al. 2001a; 
Collinge et al 2001). When described as a thermal 
component, the soft excess emission has a remark- 
ably constant temperature (~ 100 — 200 eV) across 
AGN with a wide range of black hole masses (Cz- 
erny et al. 2003; Gierlinski & Done 2004; Crummy 
et al. 2006). There are a number of ideas that have 
been proposed to explain the origin of the soft ex- 
cess component. Pounds, Done & Osborne (1995) 
pointed out the similarity of the X-ray spectrum of 
NLS1 and Cyg X-l in its high state (HS) and sug- 
gested that the soft excess could be the optically 
thick emission from an accretion disk. However, 
the observed soft excess emission is too hot for a 
disk around a supermassive black hole. One possi- 
bility is that there is a cool, optically thick Comp- 
tonizing region in addition to a hot (~ 100 keV), 
optically thin (r ~ 1) region producing the high 
energy power-law emission. The cool region can 
Comptonize the disk photons and smoothly con- 
nect the disk emission up to soft X-ray energies 
(Magdziarz et al. 1998). Gierlinski & Done (2004) 
tested this model using high quality XMM-Newton 
spectra of 26 radio-quiet PG quasars. Although 
the spectra are well described by the Comptoniza- 
tion model, the temperature of the cool Comp- 
tonizing region is remarkably constant despite a 
large range in the black hole masses. They pro- 
posed that the observed soft excess is an artifact of 
heavily smeared, strong partially ionized absorp- 
tion (see also Gierlinski & Done 2006; Schurch & 
Done 2006). Another possibility is that the emis- 
sion is dominated by ionized Compton reflection. 
This can happen if the disk is clumpy at high ac- 
cretion rates, thus hiding the hard X-ray source 
among many clumps. The ionized reflection can 
explain the observed soft excess emission (Fabian 
et al. 2002). The last two models, based on rela- 
tivistically smeared absorption and ionized reflec- 
tion, relate the soft excess emission to atomic pro- 
cesses and therefore provide a natural explanation 
for its constant temperature for AGN with a wide 
range in their black hole masses. 

The three models, cool Comptonization, ion- 



ized reflection and heavily smeared absorption, 
all result in statistically equally good fits. The 
fitting criteria are not sufficient to discriminate 
among them (Sobolewska & Done 2005). Variabil- 
ity properties of the soft and hard X-ray emission 
and the relationship between them may provide 
strong constraints on the above models for the soft 
excess emission. A 35 day long ASCA monitoring 
observation of the NLS1 galaxy Akn 564 revealed a 
distinction in the variability of the soft excess and 
power-law components on a time scale of weeks, 
with the soft excess emission varying by a factor of 
~ 6 compared to a factor of ~ 4 in the power-law 
component (Turner et al. 2001b). Similarly, along 
~ 10 days ASCA observation of the NLS1 galaxy 
IRAS 13224-3809 suggested that the soft excess 
component dominates the observed variability on 
a time scale of ~ a week, but on shorter time 
scales (~ 20000 s), the power-law component ap- 
pears to dominate the variability (Dewangan et 
al. 2002). A 12 day long ASCA observation of an- 
other NLS1 galaxy Ton S180 suggested that the 
soft excess and the power-law fluxes are not well 
correlated on time scales of ~ 1000 s (Romano et 
al. 2002). These long ASCA observations sug- 
gest two separate continuum component unlike the 
heavily smeared, ionized absorption model based 
on a single continuum component. Gallo et al. 
(2004) found evidence for alternate leads and lags 
between the 0.3 - 0.8 keV and 3-10 keV band 
emission from the NLS1 galaxy IRAS 13224-3809. 
Dasgupta & Rao (2006) reported energy depen- 
dent delay between the soft and hard X-ray emis- 
sion in Mrk 110. These observations suggest dif- 
ferent emitting regions for the soft excess and the 
power-law component from NLS1 galaxies. In this 
paper, we use XMM-Newton observations of two 
NLS1 galaxies Akn 564 and Mrk 1044 to study the 
variability properties of the soft and hard X-ray 
emission and investigate the origin of the soft ex- 
cess emission. In Section 2, we outline the XMM- 
Newton observations of Akn 564 and Mrk 1044. 
We describe the temporal analysis in Section 3 and 
spectral analysis in Section 4 and discuss the re- 
sults in Section 5. 

2. Observation and Data Reduction 

Akn 564 was observed by XMM-Newton on 
2005 January 04 (ObsID: 0206400101) for an 
exposure time of 100 ks. Mrk 1044 was ob- 
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Fig. 1. — XMM-Newton light curves and hardness ratio of Akn 564 and Mrk 1044. Left: Combined EPIC-pn 
and MOS light curves of Akn 564 in the bands 0.2 — 0.3 keV and 5 — 10 keV bands and the hardness ratio. 
The bin sizes are 512 s. Vertical dotted lines mark the dips and peaks in the 0.2 — 0.3 keV and 5 — 10 keV 
light curves that show possible evidence for a delay. Right: Combined EPIC-pn and MOS light curves of 
Mrk 1044 in the 0.2 — 0.5 keV and 2 — 10 keV bands and the hardness ratio. The bin sizes are 128 s. The 
soft and hard band light curves do not appear to be strongly correlated. 



served by XMM-Newton on 2002 July 23 (ObsID: 
0112600301) for 8 ks. The EPIC-pn and MOS 



cameras (|Striider et alj|200lt [Turner et alj|2001ah 



were operated in the small window mode using 
the medium filter during both the observations. 
The data were reduced using version 6.5.0 of the 
SAS software. Examination of the background 
rate above 10 keV showed that the observation of 
Akn 564 was partly affected by the flaring particle 
background before an elapsed time of 7.5 ks and 
this early period was therefore excluded to obtain 
a continuous exposure with steady background. 
For temporal analysis, we used all the pn and 
MOS events with patterns < 12 and a continuous 
exposure of 91.3 ks during which both the pn and 
MOS cameras operated simultaneously. There 
were three short duration (< 1000 s), small am- 
plitude background flares (maximum count rate 
~ 2 counts s -1 above 10 keV) during the observa- 
tion of Mrk 1044. The varying background contri- 
bution to the source was corrected by subtracting 
the background rate. For spectral analysis, we 
chose a count rate cut-off criterion to exclude the 
high particle background. In case of Akn 564, the 



MOS data are affected with photon pile-up. In or- 
der to achieve the best possible spectral resolution 
and to have the best spectral calibration, we chose 
the EPIC-pn data and used the good X-ray events 
(FLAG=0) with pattern zero (single pixel events) 
for Akn 564. For Mrk 1044, we used the good 
EPIC-pn events with pattern — 4 (single and 
double pixel events). We used circular regions 
with radii of 40", centered at the peak position 
of Akn 564 and Mrk 1044, to extract the source 
spectra. We also extracted background spectra 
from appropriate nearby circular regions, free of 
sources. 

We also processed the RGS data obtained from 
the long observation of Akn 564. We extracted 
the first order source and the background spectra 
by making the spatial, order selections on the pro- 
cessed event files, and calibrated by applying the 
most recent calibration data. We also used a tem- 
poral filter to exclude the high particle background 
based on a count rate cut-off criterion. We have 
used the ISIS (version 1.4.7) spectral fitting envi- 
ronment(Houck & Denicola 2000; Houck 2002) for 
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power and energy spectral analysis. Unless oth- 
erwise stated, all errors are quoted at 90% confi- 
dence level. 
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Fig. 2. — Cross correlation function between soft 
and hard X-ray light curves of Akn 564: (a) 0.2 — 
2 keV - 2 - 10 keV, (b) 0.2 - 1 keV Vs 3 - 10 keV, 
(c) 0.2 - 0.5 Vs 4 - 10 keV and (d) 0.2 - 0.3 keV 
Vs 5 — 10 keV. The time delays refer to the soft 
band. Positive peaks imply that the variations in 
the hard band lag that in the soft band. 

3. Temporal analysis 
3.1. Akn 564 

For temporal analysis, we combined the pn and 
MOS data to increase the signal-to-noise ratio. In 
order to study cross-correlation between different 
energy bands, we extracted light curves in several 
energy bands 0.2 - 0.3, 0.2 - 0.5, 0.2 - 1, 0.2 - 2, 
2-10, 3-10, 4-10 and 5-10 keV with time bins of 
64 s. Figure [T] shows the light curves of Akn 564 
in the 0.2 — 0.3 and 5 — 10 keV bands rebinned 
to 512 s. X-ray emission from Akn 564 is highly 
variable. The 0.2 — 0.3 keV band light curve shows 
a trough-to-peak variations by a factor of ~ 3.3, 
while the hard band (5 — 10 keV) shows a trough- 
to-peak variation by a factor of ~ 3.6. Although 
the 0.2 - 0.3 keV and 5-10 keV band light curves 
appear to be correlated, the most rapid variability 
is seen in the hard band. A visual comparison 
of the light curves suggests that variations in the 
hard band appear to lag those in the soft band. 



To investigate further, we have performed cross- 
correlation and power spectral analysis. 

3.1.1. Cross Correlation Analysis 

We have computed cross-correlation function 
(CCF) using the XRONOS program crosscor. 
We used the slow direct Fourier algorithm to com- 
pute the CCF between the soft and hard band 
light curves with 64 s bins. Figure [2] shows the 
CCFs as a function of time delay between the en- 
ergy bands (a) 0.2 - 2 keV and 2-10 keV, (b) 
0.2 - 1 keV and 3-10 keV, (c) 0.2 - 0.5 keV and 
4-10 keV and (d) 0.2 - 0.3 keV and 5-10 keV. 
The 0.2 - 2 keV and 2-10 keV band light curves 
are strongly correlated without any time delay as 
indicated by the sharp peak at zero lag. How- 
ever, the strength of the sharp peak decreases with 
the increase in the separation between the soft 
and hard bands. In the CCF of 0.2 — 1 keV and 
3 — 10 keV band light curves, the sharp peak weak- 
ens and another broad peak at positive time lags 
becomes relatively stronger. The broad peak dom- 
inates the CCF for the 0.2-0.5 keV and 4-10 keV 
band light curves and only a weaker sharp peak 
is seen. The sharp peak is not clearly seen in 
the CCF between 0.2 - 0.3 keV and 5-10 keV 
bands, only the broad peak is seen. Thus the hard 
4—10 keV or 5—10 keV band clearly lags the soft 
0.2 - 0.5 keV or 0.2 - 0.3 band. 

It is hard to assess the statistical significance of 
the time lag based on model fitting to the CCF 
using the traditional x 2 minimization technique. 
The reasons are (?) the shape of the CCF is gen- 
erally complex, (ii) the CCF data points are not 
independent and (Hi) the errors are not normally 
distributed. Therefore, we have adopted a model 
independent way, following Peterson et al. (1998), 
to measure the time lag. The principal source 
of uncertainty in the time lag between two long, 
continuous light curves is the uncertainty in the 
observed count rates. To account for this uncer- 
tainty, we have performed Monte Carlo simula- 
tions using IDL (version 6.3). We simulated 1000 
pairs of soft and hard band lightcurves by modi- 
fying the observed count rates by random Gaus- 
sian deviates assuming that the errors on the ob- 
served count rates are normally distributed. We 
calculated the time lag between each pair of the 
simulated soft and hard band light curves as the 
centroid of the CCF based on all points with cor- 
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Fig. 3. — Distributions of time lags obtained from the simulated soft and hard band lightcurves. Left: Time- 
lag distribution between 0.2 — 0.5 and 4—10 keV bands. Right: Time-lag distribution between 0.2 — 0.3 and 
5-10 keV bands. 



relation coefficients in excess of 0.4. Figure [3] 
shows the distribution of time lags between the 
0.2 - 0.5 keV and 4-10 keV bands. We calcu- 
lated the time lag between the observed soft and 
hard bands as the centroid of the CCF and the 
errors on the time lag by directly integrating the 
time-lag distribution over a range to obtain 68%, 
90% or 99% confidence level. Thus we measured 
the time lag between the 0.2 — 0.5 and 4—10 keV 
bands to be 1768 ± 92 s (68%), ±122 s (90%) or 
±191 s (99%). A similar procedure resulted in the 
time lag between the 0.2 — 0.3 and 5 — 10 keV to 
be 1767 ± 196 s (68%), ±320 s (90%) or ±462 s 
(99%). 

We have also searched for variations in the CCF 
or time delay. For this purpose, we divided both 
0.2 — 0.5 keV and 4—10 keV light curves into nine 
equal segments with lengths 10.1 ks and computed 
the CCFs for each segment. The nine CCFs are 
plotted in Figure [U The CCFs are not similar. 
There is no or weak correlation between the soft 
and hard band light curves in the segments II- VI 
and VIII. In segment IV, the soft and hard band 
appear to be weakly anti-correlated with a soft 
band lead of ~ 1200 s. The two bands are strongly 
correlated without any lag for the segments I & IX, 
while for the segments VII and VIII, the hard band 
lags behind the soft band. For these segments, 
the time lag varies from ~ 900 — 2900 s. In seg- 
ment VII, the soft and hard band appear to anti- 
correlate at a lag of ~ —3000 s. Brinkmann, Pa- 
padakis & Raeth (2007) have employed the sliding 
window technique to study the cross-correlation 
properties of Akn 564 between the 0.3 — 1 keV 



and 3—10 keV bands using the same data. There 
results are similar to the presented here if we note 
that the first 7.5 ks of data were excluded from 
our analysis. The leads and lags in the X-ray 
emission from Akn 564 are also similar to the al- 
ternate leads and lags in 4500 s intervals found 
between the 0.3 - 0.8 keV and 3-10 keV band 
X-ray emission from IRAS 13224-3809 by Gallo et 
al. (2004). 
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Fig. 4. — Cross correlation functions of 0.2 — 
0.5 keV and 4—10 keV band lightcurves for 
Akn 564. The CCFs have been derived by split- 
ting both 0.2 - 0.5 keV and 4-10 keV lightcurves 
into nine segments. Each segment has a length of 
10.1 ks. 
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Fig. 5. — Left: Power density spectra of Akn 564 in the full band 0.2 — 10 keV. The contribution due to the 
Poisson noise has not been subtracted. Right: Power density spectra in the soft (0.2 — 0.5 keV) and hard 
(3 — 10 keV) bands, derived from the combined EPIC-pn and MOS data. The constant power expected from 
the Poisson errors have been subtracted. The best-fitting broken power-law models are also shown. The soft 
band power density spectrum is clearly steeper and weaker than the hard band power density spectrum. 



3.1.2. The power density spectra 

To make a further study of the variability prop- 
erties of Akn 564, we have calculated power den- 
sity spectra (PDS) of Akn 564 from the combined 
pn+MOS lightcurves in the full (0.2-10 keV), soft 
(0.2 - 0.5 keV) and hard (3 - 10 keV) bands, sam- 
pled at 16 s. We used the slow algorithm imple- 
mented in the XRONOS program 'powspec'. Fol- 
lowing Papadakis & Lawrence (1993), we rebinned 
the PDS in logarithmic space with a binsize of 20 
and performed the fitting within the ISIS (version 
1.4.5) spectral fitting environment. We converted 
the PDS to equivalent energy spectra (power — > 
counts/bin; Hz — > keV) using the sitar timing 
package. We fitted the PDS with a simple power- 
law model. We also used a constant to account for 
the power arising from the Poisson noise. The sim- 
ple power-law model resulted in an unacceptable 
fit (minimum x 2 = 225.2 for 139 degrees of free- 
dom (dof)). We then replaced the power-law with 
a broken power-law model. This model resulted in 
the minimum x 2 = 169.2 for 137 dof. The broken 
power-law model provided an improved fit at a sta- 
tistical significance level of > 99.99% compared to 
the power-law model according to the maximum 
likelihood ratio test. Thus we conclude that the 
PDS of Akn 564 has a break at (1.2±Q.3)xlO" 3 Hz 
with power-law indices I.S^q'J, below the break 



and 3.2^0 4 above the break. The break frequency 
reported here is similar to the high frequency 
break at ~ 2 x 10~ 3 Hz discovered by Papadakis 
et al. (2002) based on ASCA data (see also Ale- 
varo et al. 2006) The constant power due to Pois- 
son noise is 0.024±0.0012 (rms/mean) 2 Hz- 1 . The 
PDS and the best-fit broken power-law plus con- 
stant model are plotted in Figure [5J The power 
arising from the Poisson errors have not been sub- 
tracted. 

We also fitted the power-law and broken power- 
law models to the soft (0.2—0.5 keV) and hard (3— 
10 keV) PDS of Akn 564. The power-law model 
resulted in an unacceptable fit to the PDS in the 
soft band (x 2 = 193.2 for 139 dof) with power- 
law index of 1.7 ± 0.1. The broken power-law 
model provided a better fit (x 2 /dof — 160.9/137) 
with power-law indices 1.4^' 2 and 4.0l^'g before 
and after the break frequency l.lj^'g x 10 -3 Hz. 
This is an improvement over the power-law fit 
at a significance level of > 99.99% based on the 
maximum likelihood ratio test. The power-law 
model provided an acceptable fit to the hard band 
PDS (x 2 = 138.4 for 139 dof) with an index of 
1.2 ± 0.1. The broken power-law model resulted in 
an improved fit (x 2 /dof = 127.2/137) with indices 
< 1.0 (90% upper limit) and 1.8±g | before and 
after the break frequency < 1.2 x 10~ 3 Hz (90% 
upper limit). This is an improvement over the 
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power-law model at a level of 99.6% significance 
level. However, the parameters are not well deter- 
mined. We have plotted the soft and hard band 
PDS and the best-fit broken power-law models in 
Fig. [5] (right panel). Clearly the soft band PDS 
is steeper than the hard band PDS. This implies 
that the X-ray emission in the hard band varies 
more than that in the soft band on shorter time 
scales. 
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Fig. 6. — Cross correlation function between soft 
and hard X-ray light curves of Mrk 1044: (a) 0.2 — 
2 keV - 2 - 10 keV, (b) 0.2 - 1 keV Vs 3 - 10 keV, 
(c) 0.2 - 0.5 Vs 4 - 10 keV and (d) 0.2 - 0.3 keV 
Vs 5 — 10 keV. The time delays refer to the soft 
band. Positive peaks imply that the variations in 
the hard band lag that in the soft band. 

3.2. Mrk 1044 

As stated above, the observation of Mrk 1044 
was affected by short-duration, weak flares. There- 
fore, the source light curves were corrected for 
background contribution using appropriate back- 
ground light curves extracted from blank sky re- 
gions. The combined pn+MOS and background 
corrected light curves of Mrk 1044 are shown in 
Fig. Q] (right). X-ray emission from Mrk 1044 
is variable with trough-to-peak variations of ~ 2 
(0.2 - 0.5 keV) and - 3 (4 - 10 keV). However, 
the 0.2 - 0.5 keV and 4-10 keV light curves do 
not appear to be strongly correlated. 

The cross correlations between the soft and 
hard bands are plotted in Figure El Similar to 



Akn 564, the 0.2 - 2 keV and 2-10 keV bands 
are strongly correlated without any delay and the 
strength of this correlation decreases with the in- 
terval between the soft and hard bands. There 
is no correlation between the 0.2 — 0.3 keV and 
5 — 10 keV bands at zero lag. There is a weak 
anti-correlation at a positive delay of ~ 1000 s 
and a weak correlation at positive delays of 3000 — 
5000 s. As the CCF for Mrk 1044 is based on a 
short 7.9 ks observation, the CCF is likely only one 
representation of many possible cross correlations 
similar to that seen for Akn 564. 

4. Spectral Analysis 

We have performed spectral modeling of Akn 564 
and Mrk 1044. We present the results based on 
the high signal-to-noise pn data for both Akn 564 
and Mrk 1044 and the high resolution RGS data 
for Akn 564. We note that Papadakis et al. (2007) 
have performed a detailed spectral analysis of 
Akn 564 using the same observation. The pur- 
pose here is to investigate the alternative models 
for the origin of the soft X-ray excess emission. 
All the spectra were analyzed with the Interactive 
Spectral Interpretation System (ISIS, version 
1.4.5). The pn data were grouped to a mini- 
mum of 100 and 20 counts per spectral channel 
for Akn 564 and Mrk 1044, respectively. We con- 
sider a model fit as statistically unacceptable if 
the the null hypothesis can be rejected at more 
than 95% significance. The errors on the best-fit 
spectral parameters are quoted at a 90% confi- 
dence level. Below we use narrow Gaussian line 
components with their widths fixed at a = 10 eV 
that are unresolved by the EPIC-pn data. 

4.1. The soft X-ray excess emission 

NLS1 galaxies are known to show strong soft X- 
ray excess emission below ~ 2 keV. To show this 
soft excess, first we fitted a simple absorbed power 
law (PL) model to the EPIC-pn spectra above 
3 keV (3 - 10 keV for Mrk 1044 and 3 - 11 keV 
for Akn 564). The simple PL model resulted in 
minimum x 2 = 527.4 for 519 dof and 119.0 for 
127 dof for Akn 564 and Mrk 1044, respectively, 
thus providing good fits to the 3 — 10 keV pn data. 
We have plotted the ratio of the pn data and the 
best-fit PL model in Figure [7] after extrapolating 
the model to low energies. The plots clearly show 
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strong soft excess emission below ~ 2.5 keV and 
possible iron Ka lines near ~ 6.4 keV in the spec- 
tra of both the NLS1. The soft X-ray excess emis- 
sion from NLS1 are usually well described by a 
simple blackbody or a multicolor disk blackbody 
(MCD). Addition of the MCD component to the 
absorbed PL improved the fit. The fit was reason- 
ably good for Mrk 1044 (x 2 /dof = 628.0/541). 
Addition of a narrow Gaussian line at ~ 6.4 keV 
further improved the fit {x 2 /dof = 611.7/539). 
The best-fit parameters of the MCD+PL model 
for Mrk 1044 are listed in Tableland the EPIC-pn 
data, the best-fit model and residuals are plotted 
in Figure [9l 
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Fig. 7. — Ratios of observed EPIC-pn data for 
Akn 564 and Mrk 1044 and the corresponding 
best-fit power-law model fitted above 3 keV and 
extrapolated to lower energies. 

The MCD+PL model was a poor fit for Akn 564 
(x 2 /dof = 1494.1/913). Examination of the resid- 
uals showed broad emission/absorption features 
below ~ 1 keV in the spectrum of Akn 564. To 
verify and investigate these features, we have an- 
alyzed the high resolution soft X-ray spectra of 
Akn 564 obtained with the two RGS simultane- 
ously with the EPIC-pn data. We used the Cash 
statistic to fit the low count RGS data. We fit- 
ted a simple PL model, modified with the neu- 
tral absorption, jointly to RGS1 and RGS2 data 
in the 0.35 — 1.9 keV band. This resulted in a steep 
spectrum (r ~ 3.0) and N H = 7.6 x 10 20 cm" 2 . 
Addition of an MCD component improved the 
fit (C/dof = 2575.9/1572; AC* = -225.4 for 
two parameters). The best-fit parameters are 
N H ~ 5.2 x 10 20 cm" 2 , T ~ 2.44 and the MCD 
kTi n ~ 165 eV. Examination of the ratio of the 



data and model revealed a strong absorption fea- 
ture at ~ 0.5 keV and additional weak and narrow 
absorption features in the 0.5 — 1 keV band, sug- 
gesting the presence of a warm absorber medium 
in Akn 564. 

4.2. Warm absorber model for Akn 564 

To describe the weak absorption features seen 
in the RGS data, we have created a warm absorber 
model for Akn 564 using the spectral simulation 
code CLOUDY version 7.02.01 (last described by 
Ferland et al. 1998). We used a multicomponent 
ionizing continuum similar to that observed from 
Akn 564. We used the 'agn' continuum available 
within CLOUDY. This continuum has four param- 
eters: the temperature of the big blue bump emis- 
sion, the slope of the high energy X-ray contin- 
uum, X-ray to UV ratio (a ox ) and the UV spec- 
tral index. We set the temperature of the big blue 
bump component and the slope of the high en- 
ergy continuum to be similar to the temperature 
of the soft X-ray excess emission and the slope 
of the hard power law, respectively, derived from 
the EPIC-pn data. For Akn 564, a ox = 0.941 
(Gallo 2006) and we used a typical value of —0.5 
for the UV spectral index (Elvis et al. 1994; Fran- 
cis 1993). The X-ray power law is assumed to fall 
as v z above 100 keV and the big blue bump com- 
ponent is assumed to have an infra-red exponential 
cut-off at kT IR = 0.01 Ryd or 0.136 eV. We also 
included the cosmic microwave background radia- 
tion so that the incident continuum has nonzero 
intensity at very long wavelengths. We assumed 
a plane parallel geometry and calculated grids of 
models by varying the ionization parameter and 
the total hydrogen column density. We also in- 
cluded UTA features from Gu et al. (2006) in our 
calculation. The grid of models were imported to 
ISIS in the form of an XSPEC-stylc multiplicative 
table model as described in Porter et al. (2006). 

To describe the weak absorption features in 
the RGS spectra of Akn 564, we multiplied the 
CLOUDY warm absorber model to the MCD+PL 
model. The model was also modified by neu- 
tral absorption along the line of sight as be- 
fore. The addition of the warm absorber com- 
ponent improved the fit (C/dof = 2469.3/1569; 
AC = —106.6 for three additional parameters) 
with warm absorber column, N\y ~ 4xl0 20 cm -2 , 
ionization parameter {log^ ~ 2) and the out- 
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Table 1: Spectral model parameters for Akn 564 and Mrk 1044. 



Parameter 




Akn 564 1 


Mrk 1044 


MCD+PL 


Nh (10 2U cm" 




5.24(f) 


3.54(f) 


kT in (eV) 




179 ±2 


138 ±3 


TIMCD 

r 




300211°* 
2.48 ± 0.01 


2992± 3 ^ 
2.14 ±0.04 


riPL 3 




(1.35 ±0.02) x 10~ 2 


(3.2 ±0.1) x 10~ 3 


E FeKa (keV) 
f 4 

JFeKa 

f 5 
J0.3-10keV 




6.72 ± 0.03 


6.42 ± 0.04 




(5.1 ± 1.8) x 10~ 6 
6.6 x 10~ n 


(1.14 ±0.46) x 10~ 5 
2.1 x 10~ n 


X 2 /dof 




919.3/854 


611.7/539 


Blurred Reflection 


N H (10 2U cm" 

r 

riPL 3 
P 




5.24(f) 
2.4681 ,;™ 2 
9.29 ±0.03 x 10~ 3 
> 10.0 


3.54(f) 

9 9n +0.002 
z - zu -0.004 

1.43±g;g| x 10~ 3 
> 10.0 


i 




< 1.58 

66 ± 1° 


< 1.25 

65.7«; 3 ° 


Tlrefl 6 

£ (erg cm s _1 
E FeKa (keV) 
f 4 

JFeKa 

f 5 


) 


(4.27 ±0.02) x 10~ 7 
0.69 ±0.01 
34521 1 " 

n 79 +0.05 
D -' z 0.03 

2.5t 2 '.5 x fO -6 
6.4 x 10" 11 


1.7+° ° 2 x 10" 6 
H5+ 04 
334.5l 2 ;2 
6.42 ± 0.04 

i.olg:S x io~ 5 

2.0 x 10" 11 






954.2/850 


671.6/535 


Smeared absorption 


N H (10 2U cm" 

r 




5.24(f) 

9 7QO+0.006 

Z. IOd_Q QQg 


3.54 (f) 
2.49 ±0.02 


npz, 3 

wind JVh (10 22 cm" 2 ) 
log (£/ erg cm s _1 ) 
a 7 


(3.67 ±0.05) x 10~ 2 

9 Q+0 05 
z ' a -0.03 

0.81 ±0.03 


(7.6 ±0.2) x 10~ 3 

2.89 ±0.04 
0.29 ±0.03 


E F eKa (keV) 
JFeKa 




6.72 ± 0.03 


6.42 ± 0.04 




(5.1 ±2.6) x 10~ 6 


(1.2 ±0.5) x 10~ 5 


f 5 
/0.3-lOfceV 




7.1 x 10~ n 


2.1 x 10~ n 






1147.1/853 


592.4/538 


Thermal Comptonization 


N H (10 2U cm" 

r 




5-24(/) 
2.46 ±0.01 


3.54 (f) 
2.11 ±0.05 


riPL 3 




(1.31 ±0.02) x 10~ 2 


(3.1 ±0.15) x 10~ 3 


kT seed { eV) 
kT e (keV) 




30(f) 
0.176 ±0.006 


30(f) 
n 1 07+0.012 


r thcomp 




2.1 ± 0.1 


2.1 ± 0.2 


E F eKa (keV) 

f 4 
JFeKa 




6.72 ± 0.03 


6.42 ± 0.04 




(4.9 ± 1.8) x 10~ 6 


(1.1 ±0.5) x 10~ 5 


t 5 
/0.3-lOfceV 




6.7 x 10~ n 


2.1 x 10~ n 






900.2/853 


603.0/538 



1 For Akn 564, all models were modified by the two phase warm absorber medium inferred from the RGS data and the EPIC-pn 

data were fitted in the 0.6 — 11 keV band. 
2 MCD normalization n^CD = (Rin / km) / (D / 10 kpc), where Ri n is the inner radius and D is the distance. 
3 Powcr-law normalization in units of photons keV -1 cm~ 2 s _1 at 1 keV. 
4 ine flux in photons keV -1 cm -2 s _1 . 

5 Obscrvcd flux in units of ergs cm~ 2 s — 1 in the 0.3 — 10 keV band. 
6 Normalization of the reflected spectrum. g 
7 Gaussian sigma for velocity smearing in units of v/c. 
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Fig. 8. — (ab) RGS1 data, the best-fitting joint model and the ratios of the observed data and the model 
for Akn 564. (b) The same as (a) but for RGS2. The best-fit model was derived by joint spectral fitting to 
both RGS1 and RGS2 data. 



flow velocity (v ~ 300 km s _1 ). An additional 
warm absorber component further improved the 
fit (AC = —70.5 for three parameters) with 
N w ~ 2 x 10~ 20 cm" 2 , logS, < 0.3 and v ~ 



-1000 km s" 



Examination of the residuals 



showed a strong absorption feature near 0.5 keV 
at the location of the neutral oxygen edge. How- 
ever, the shape of the absorption feature does not 
resemble to an absorption edge, it is similar to an 
absorption line with a very weak emission line at 
the center. This feature is likely due to possible 
RGS calibration errors at the oxygen edge and/or 
from the Galactic oxygen along the line of sight. 
Addition of a Gaussian absorption line improved 
the fit (AC = —97.3 for three parameters). The 
parameters of the best-fit model to the RGS data 
is listed in Table [2l We note that the best-fit neu- 
tral absorption column density derived from the 
RGS data is N H = 5.24^°;^ x 10 20 cm" 2 which is 
less than the Galactic column of 6.40 x 10 20 cm -2 
derived from the HI map of Dickey & Lockman 
(1990). The best-fit column, however, is consis- 
tent with the Galactic column of 5.34 x 10 20 cm -2 
derived from the Leiden/ Argentine/Bonn (LAB) 
Survey of Galactic HI (Kalbera et al. 2005). McK- 



ernan, Yaqoob, &; Reynolds (2007) derived cold 
Nh = 4^ 2 x 10 20 cm -2 using the Chandra high en- 
ergy transmission grating observations of Akn 564. 
Given the above uncertainties in the measurement 
of cold Nh, we fix the cold absorption at the best- 
fit RGS value of Nh = 5.24 x 10 20 cm" 2 in all the 
subsequent fits for Akn 564. 

The above analysis clearly shows the presence 
of a two phase warm absorber medium in Akn 564. 
Some or all of the emission/absorption features re- 
sulting in the poor quality of the fit to the EPIC- 
pn, described above in section 4.1 are likely due 
to the presence of warm absorbers. However, it is 
difficult to determine accurately the parameters of 
the warm absorber based on the lower resolution 
EPIC-pn data alone. Therefore, we used the best- 
fit warm absorber components derived from the 
RGS data to model the EPIC-pn data. We used 
the two component warm absorber model with the 
fixed parameters along with the MCD+PL model 
and fitted the EPIC-pn data. The addition of the 
warm absorber components improved the fit, pro- 
viding x 2 /dof — 1654.1/912. However, the fit is 
still not acceptable mainly due to a strong ab- 
sorption feature near 0.5 keV and slight excess of 
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Table 2: Best- fit model parameters derived from the high resolution RGS observation of Akn 564 
Neutral absorption N H (10 2U cm" 2 ) 5.25^;^ 

Warm absorber I N w (10 20 cm- 2 ) 3.9±g:| 
Zog£ 2.0 ± 0.1 

(u/kms- 1 ) 1 — 280±i|g 



Warm absorber II N w (10 20 cm" 2 ) 2.1 



+0.2 
-0.1 



logi < 0.32 

1 

-155 



(v/kms- 1 )* -1080t 175 



Disk blackbody kT in (eV) 171 ± 2 

n MCD b 29171^ 

Power law T 2.48 ± 0.01 

n Pi c (1.50 ± 0.01) x 10" 2 

Absorption line E (keV) 0.541 ± 0.001 

fune d ' -2-61^ X 10- 4 

C/dof 2301.5/1577 



"Velocity with respect to the systemic velocity. Negative sign indicates outflow. 

'MCD normalization umcd = {Rin I 'km) / (D / '10 kpc), where Ri n is the inner radius and D is the distance. 
c Power-law normalization in units of photons keV - 1 cm -2 s _1 at 1 keV. 
d ine flux in photons keV~ 1 cm -2 s _1 . 



emission below 0.5 keV. The absorption feature 
is stronger by a factor of ~ 2.5 in the EPIC-pn 
data than that derived from the RGS data. Also 
the excess emission below 0.5 keV is not evident 
in the RGS data. Thus EPIC-pn and RGS do 
not agree well below 0.6 keV. Therefore we ignore 
the EPIC-pn data below 0.6 keV in the fits for 
Akn 564 described below. 

After excluding the data below 0.6 keV, the 
MCD+PL model modified with the two com- 
ponent warm absorber and neutral absorption 
along the line of sight provided a much better fit 
(x 2 /dof = 941.2/856). Adding a narrow Gaus- 
sian line at ~ 6.6 keV further improved the fit 
(Achi 2 /dof — —21.9 for two parameters). Fig- 
ure [TO] shows the EPIC-pn data, the best-fit model 
and the x-residuals and the best-fit parameters are 
listed in Table Q] The parameters of the MCD and 
PL components are similar to that derived from 
the RGS data (see Table [2 and |2J). The MCD com- 
ponent used above to describe the soft X-ray ex- 
cess emission is physically inappropriate. The in- 
ner accretion disk temperatures, kTi n — 179±2 eV 
for Akn 564 and 138 ±3 eV for Mrk 1044 are much 



higher than that expected from a disk around a 
black hole of mass 10 6 — 10 7 M Q . For this reason, 
the origin of the soft excess emission has remained 
unsolved for more than two decades. 

4.3. Relativistically blurred ionized reflec- 
tion 

One possibility is that the reflected emission 
from a partially ionized accretion disk could be 
strong below 1 keV and can account for the strong 
soft excess emission (see e.g., Crummy et al. 
2005). Though there is no clear evidence for a 
broad iron Ka line, such a line from an accre- 
tion disk could be hidden in the strong X-ray 
continuum. To test whether the strong soft X- 
ray excess emission from Akn 564 and Mrk 1044 
could be produced by the ionized reflection emis- 
sion arising from the reprocessing of the primary 
X-ray emission in the disk, we have fitted the rela- 
tivistically blurred ionized reflection model to the 
EPIC-pn spectra of both the NLS1. We used the 
table model ref lion (Ross & Fabian 2005; Ross, 
Fabian & Young 1999) for the ionized reflection 
and blur the ionized reflection emission relativisti- 



11 




cally by convolving with a LAOR kernel (kdblur) to 
obtain the relativistically broaden ionized reflec- 
tion from the inner regions of an accretion disk. 
The parameters of ref lion model are ionization 
parameter (£), iron abundance relative to solar 
(Ap e ), photon index of the illuminating power- 
law and the normalization of reflected spectrum. 
The parameters of the convolution model kdblur 
are the inner and outer radii of the disk (r» n and 
r out ), emissivity index {(3) and the disk inclina- 
tion (i). We fixed the outer radius at r out — 400r g 
where r g = GM/c 2 is the gravitational radius. 
We used the power-law component for the con- 
tinuum. We also used a narrow Gaussian to de- 
scribe the unresolved iron Ka line from distant 
matter. We also used the two component warm 
absorber model with the parameters fixed at the 
values derived from the RGS data. We refer the 
full model as the blurred reflection model. The 
best-fit parameters of this model are listed in Ta- 
ble[JJ The observed data and the best-fit reflection 
models are shown in Figure [9] (Mrk 1044) and Fig- 
ure [10] ( Akn 564) . The model resulted in poorer 
fit than the MCD+PL model for both the NLS1 



(x 2 /dof = 954.2/837 for Akn564 and 671.6/535 
for Mrk 1044). Also the emissivity index is un- 
physically large (/3 > 10) and the inner radius is 
very small for both the AGN, suggesting that the 
model cannot account for the observed smooth- 
ness of the soft excess X-ray emission. The very 
steep emissivity indices and the small inner radii 
of the LAOR kernel were required to smooth the 
recombination features sufficiently to produce the 
featureless soft X-ray excess emission from both 
the NLSls. 

4.4. Smeared absorption from partially 
ionized material 

The observed soft excess emission could be an 
artifact of smeared absorption from partially ion- 
ized material as described by Gierlin'ski & Done 
(2004). To test this scenorio, we have fitted the 
EPIC-pn data with the smeared absorption model. 
We used the XSPEC local model swindl for the 
smeared absorption from a partially ionized wind 
from an accretion disk and the simple PL model 
for the continuum. Again we used a Gaussian 
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Fig. 10. — Observed EPIC-pn data and the best-fitting MCD+PL, blurred ionized reflection, smeared 
absorption from partially ionized wind and optically thick thermal Comptonization models for Akn 564. 
The best-fit parameters are listed in Table [TJ 



for the narrow iron Ka line and the two compo- 
nent warm absorber model with fixed parameters 
derived from the RGS data. The parameters of 
the smeared wind model are absorption column 
density, the ionization parameter (£ = L/nr 2 ) 
and the Gaussian sigma for velocity smearing in 
units of v/c. The smeared absorption model pro- 
vided a good fit to the spectrum of Mrk 1044 
(x 2 /dof = 592.4/538) but fails to describe the 
spectrum of Akn 564 (x 2 /dof = 1147.1/842). 
This model requires unphysically large smearing 
velocity (~ 0.8c) for Akn 564 to smooth suffi- 
ciently the absorption lines in order to produce 
the observed smooth soft excess emission (see Ta- 
ble OQ and Fig.HDJ). 

4.5. Optically-thick, thermal Comptoniza- 
tion 

The strong and smooth soft excess emission 
from NLS1 can also be described as arising from 
thermal Comptonization in an optically thick 
corona. To test such a scenario, we used the ther- 
mal Comptonization model (ThComp) described 



by Zdziarski et al (1996). The free parameters 
of ThComp model are the asymptotic power-law 
index {Tthcomp), electron temperature (kT e ) and 
the seed photon temperature {kT^). The electron 
scattering optical depth (r) can be calculated from 
the asymptotic power-law photon index (F compt h) 
and electron temperature (kT e ) as follows 



1 



"I 1/2 



(kT e /m e c 2 )T(l+T/3) 



- (1) 

2 V ; 



with Tthcomp — a + 1 (Sunyaev & Titarchuk 1980). 
We also used the PL model to describe the hard 
PL for the hot Comptonized component. We fixed 
the temperature of the seed photons for the op- 
tically thick Comptonization at 30 eV. This is 
equivalent to assuming that the seed photons for 
the Comptonizations in the warm and hot plasma 
arise from the accretion disk. Again we used 
the narrow Gaussian line for the iron Ka line. 
This model provided good fits for both the NLSls 
(x 2 /dof = 900.2/842 for Akn 564 and x 2 /dof = 
603.0/538 for Mrk 1044). The best-fit parameters 
are listed in Table Q] and the EPIC-pn data and 
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the best- fit model are plotted in Fig. [9] and Fig.flOl 
. The electron temperatures are low (~ 0.18 keV 
for Akn 564 and - 0.14 keV for Mrk 1044) and 
electron scattering optical depths are very high 
(t - 45) for both the NLSls. 

We note that the thermal Comptonization 
model best describes the EPIC-pn data in the 
0.6-11 keV band for Akn 564 (see Table CQ). 
We expect that the most suitable model must 
provide a reasonably good fit to the full band 
EPIC-pn data except for the absorption feature 
at 0.54 keV, provided that the calibration errors 
at low energies (< 0.5 keV) are at a level of a 
few percent. We find that the thermal Comp- 
tonization model best describes the full band 
data with x 2 /do/ = 1164.8/911 without any sig- 
nificant change in the best-fit parameters listed 
in Table Q] Addition of an absorption line at 
0.54 keV improved the fit to x 2 /dof = 985.5/908 
with line center at E — 0.54 ± 0.005, fu ne 

: -",,1-1 ± 0.06) x 10~ 4 ergs cm" 2 s" 1 and width 
cr < 8 eV. The ionized reflection and smeared 
absorption models both fail to describe the data 
below 0.5 keV when extended to lower energies. 
Refitting the reflection model with the addition 
of an absorption line at 0.54 keV to the full band 
data resulted in poor fit {x 2 /dof = 1392.2/906) 
with large emissivity index (/? = 9.17igQg) and 
small inner radius (r^„ < 1.239r s ) as before. The 
smeared absorption model with the addition of 
the absorption line provided the worst fit among 
the four models. 

5. Discussion 

We have studied the X-ray temporal and spec- 
tral characteristics of Akn 564 and Mrk 1044. We 
found that the 0.2 - 2 keV and 2-10 keV band 
lightcurves are strongly correlated without any 
time delay. However, the strength of this correla- 
tion decreases with increasing separation between 
the soft and hard bands. We have found signifi- 
cant correlation between the soft (0.2 — 0.3 keV, 
0.2 - 0.5 keV) and hard (5 - 10 keV, 4-10 keV) 
band lightcurves at positive time delays. The 
4-10 keV band lags behind the 0.2 - 0.5 keV 
band by 1767 ± 122 s (at 90%) in Akn 564. We 
note that Arevalo et al (2006) found time delays 
that increase both with time scale and energy sep- 
aration between the energy bands. This means 



that there is no single, well defined delay. Instead, 
the measured delay based on the cross-correlation 
analysis is the averaged delay between the soft and 
hard bands. We find no correlation at zero time 
lag between the 0.2 — 0.3 keV and 5 — 10 keV band 
emission from Mrk 1044. We confirm the pres- 
ence of a break at ~ 1.2 x 10~ 3 Hz in the power 
density spectrum of Akn 564, earlier detected by 
Papadakis et al. (2002) based on ASCA observa- 
tions. We also find that the PDS of Akn 564 is en- 
ergy dependent. The hard band (3 — 10 keV) PDS, 
consistent with a simple power law, is flatter and 
stronger than the soft band (0.2 - 0.5 keV) PDS 
that shows clear evidence for the break. These ob- 
servations strongly suggest that different variabil- 
ity process dominate the soft and hard band emis- 
sion. This may suggest physically distinct spectral 
components for the soft excess and the power-law 
components. 

5.1. Implications of the observed time de- 
lays 

The observed hard band delays from Akn 564 
demonstrate that the soft excess emission cannot 
be the reprocessed emission of the primary X-ray 
emission. Fabian et al. (2002) presented the re- 
flection dominated ionized accretion disk model 
and successfully reproduced the strong soft excess 
emission from the NLS1 galaxy 1H 0707-495. In 
this model, the strong soft excess is caused by the 
multiple reflections of mostly hidden hard X-ray 
continuum by the clumpy disk and the stronger 
reflectivity of the partially ionized material be- 
low ~ 0.7 keV, associated with an abrupt change 
in the opacity (Ross, Fabian & Ballantyne 2002; 
Fabian et al. 2004). The soft X-ray excess emis- 
sion is the reprocessed emission in this model, 
hence it is unlikely that the soft X-rays lead the 
hard X-ray emission. This is, however, not a 
strong constraint on the ionized reflection model. 
Many of the soft X-ray lines expected from the 
accretion disk are excited by soft X-ray photons 
with energies slightly greater than the line ener- 
gies as the photoelectric absorption cross section 
varies as ~ E~ 3 . The ionizing soft X-ray pho- 
tons can lead the hard (3 — 10 keV) band emis- 
sion due to the Compton delay of hard photons in 
the Comptonization process. Since the life times 
of excited states, responsible for the line photons, 
are very short compared to the expected Compton 
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delay between the soft and hard photons, the line 
photons as well as the scattered soft X-ray emis- 
sion can still lead the hard band primary emission. 
However, in this case, the soft X-ray excess emis- 
sion is not expected to show rapid and large am- 
plitude variability as the reflection component and 
broad iron line arc known to show little variability 
in MCG-6-30-15 (Miniutti et al. 2007). Moreover, 
the blurred ionized reflection model with physi- 
cally plausible parameters cannot reproduce the 
observed soft X-ray excess emission from Akn 564 
and Mrk 1044. 

Gierlinski & Done (2004) presented an alterna- 
tive model for the origin of the soft excess emis- 
sion. In this model, the soft excess emission is 
mainly due to the strong and smeared absorp- 
tion near ~ 0.7 keV. The soft excess and the 
hard power law are part of the same and single 
continuum component except for the small contri- 
bution from smeared emission lines. This means 
the soft and hard band emission arise from the 
same continuum component or physical process - 
Comptonization. The time delay is a natural out- 
come of the Comptonization process. A general 
prediction of Comptonization models is that the 
hard X-ray variations should lag behind those in 
softer bands (Payne 1980). Therefore it is natu- 
ral to attribute the observed time delays between 
hard and soft photons to this process. The time 
lag represents the difference in the photon escape 
times between the soft and hard bands. The hard 
band photons tend to have undergone more scat- 
tering events, therefore travel longer and escape 
the corona later than the soft band photons. Thus 
the observed time lag is consistent with a single 
continuum component and hence with the smeared 
absorption model. However, there are three prob- 
lems with the smeared wind model, (i) This model 
(XSPEC local model swindl; Gierlinski & Done 
2004) does not describe the X-ray spectrum of 
Akn 564 satisfactorily, (ii) The 3 — 10 keV band 
shows more variability power on short time scales 
than the 0.2 — 0.5 keV band, which is not expected 
in the simple Comptonization models in which the 
variability is caused by the variations in the seed 
photons alone. (in) The smearing velocity for 
Akn 564 suggested by the swindl model is very 
large (v ~ 0.8c) for disk winds. Hydrodynamical 
simulations show that the line driven winds do not 
have the required smearing velocity (N. J. Schurch, 



private communication). It remains to be seen if 
the magnetically driven winds can provide the re- 
quired smearing. 

The two component model, Comptonization in 
two separate regions, is a viable option. In this 
model, the soft excess emission arises from the 
Comptonization in the low temperature (warm), 
optically thick plasma while the hard power-law 
is produced by the Comptonization in the hot 
plasma. Indeed the EPIC-pn spectra of Akn 564 
and Mrk 1044 are both well described by the two 
component model. The temperature of the warm 
plasma is ~ 175 eV for Akn 564 and ~ 135 eV 
for Mrk 1044. This region, giving rise the soft 
excess emission, is optically thick (t ~ 45). Phys- 
ically, this component could be the warm skin on 
the accretion disk surface or a region between the 
disk and an optically thin inner flow (Magdziarz 
et al. 1998) or an optically thick corona coupled 
with a truncated disk (Done & Kubota 2006). De- 
pending on the geometry of the disk and the two 
component corona system, there can be delay be- 
tween the soft excess and high energy power-law 
components. 

However, the two component Comptonization 
model has difficulty in explaining the constancy of 
the temperature of the Comptonizing region pro- 
ducing the soft excess emission. Gierlinski & Done 
(2004) derived the temperature of the putative 
Comptonizing region in the range 0.1 — 0.2 keV 
which is remarkably constant for 26 different PG 
quasars. They concluded that the constancy of 
the temperature is inconsistent with that expected 
from the span of the asymptotic power-law index 
(Tthcomp ~ 1-5 - 2.5). However, the 0.3 - 10 keV 
XMM-Newton band only includes the high en- 
ergy tail of the optically thick Comptonized emis- 
sion. Both the peak of the soft excess compo- 
nent (kT ~ 0.1 — 0.2 keV) and the asymptotic 
power-law component fall below 0.3 keV. Using 
such data, it may be difficult to measure accu- 
rately both the temperature and the asymptotic 
photon index, particularly in the presence of line 
of sight neutral absorption and complex warm ab- 
sorber, both modifying the soft X-ray spectrum 
below 1 keV. NLS1 galaxies show a large diver- 
sity in their soft X-ray ROSAT photon indices 
(r ~ 1.5 — 5.0), suggesting varying shape of the 
soft excess component (Boiler, Brandt & Fink 
1996). The constancy of the temperature for the 
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soft excess emission is required to be investigated 
using a sample of NLS1 galaxies similar to that 
used by Boiler, Brandt & Fink (1996). Another 
possibility is that the disk/corona geometry is dif- 
ferent in NLS1 galaxies as they are accreting close 
to the Eddington rate. In such situations, the in- 
ner disk is likely thick and hotter. It remains to 
be seen if such thick structure can have uniform 
temperatures across AGN with a large range in 
their luminosity and black hole masses. 

5.2. The disk-corona geometry 

Inverse Compton scattering by thermal elec- 
trons provides the best possible explanation for 
the observed smooth soft X-ray excess emission 
from Akn 564 and Mrk 1044. Our spectral analysis 
favors two physically distinct corona: high optical 
depth, low temperature corona and high temper- 
ature, low optical depth corona. There are two 
possible geometries with two component corona. 
The cool, optically thick corona may be similar 
to the coupled disk corona geometry proposed by 
Done & Kubota (2006) for the very high state of 
black hole binaries. In this geometry, the inner 
accretion disk is truncated and is coupled to an 
optically thick corona, and gravitational energy re- 
lease powers both the disk and the optically-thick 
corona. Comptonization of the disk photons in 
this corona may give rise to the broad soft ex- 
cess component. In addition to the optically thick 
corona, a hot corona is responsible for the power- 
law component. The hot corona could be the same 
as the patchy corona above the disk and is likely 
powered by magnetic flares in the disk similar to 
the solar flares (Haardt, Maraschi, & Ghiscllini 
1994; Stern ct al. 1995; Poutanen & Svensson 
1996; Poutanen & Fabian 1999; Merloni & Fabian 
2001). The seed photons could either be the soft 
photons from the disk or those arising from the 
cool, optically thick corona. Since the disk is cou- 
pled to the optically thick corona, appreciable time 
lag between the disk photons and the soft excess 
photons are not expected. The variability of the 
disk flux and the soft excess emission is governed 
by the variations in the accretion rate, while the 
variability of the power-law component is governed 
both by the variation in the soft photons as well 
as any variability process intrinsic to the corona 
such as that related to variation in the electron 
population. The flatter power density spectrum 



of Akn 564 in the hard band may be related to 
the additional variability of the power-law compo- 
nent due to the changes in electron population on 
short time scales e.g., caused by magnetic flares 
above the accretion disk. 

Another possible geometry is the ionized sur- 
face of a disk acting as an optically thick and cool 
corona with a temperature of ~ 0.1 — 0.3 keV (e.g, 
Magdziarz et al. 1998; see also O'Brian et al. 
2001; Page et al. 2002; Vaughan et al. 2002). In 
this case, thermal Comptonization of disk photons 
by thermal electrons in the hot surface of the disk 
will produce the smooth soft excess. The steep 
power-law extending to high energies can be pro- 
duced in another purely hot thermal/nonthermal 
or in a hybrid thermal-nonthermal plasma. 

It is difficult to distinguish the two geometries 
described above. The different power spectra of 
the soft and hard band X-ray emission may con- 
strain the two geometries. Recently McHardy et 
al. (2006) have shown a strong correlation be- 
tween the break frequency, black hole mass and the 
bolomctric luminosity for black hole binaries and 
AGN. If the break frequencies are related to the 
characteristic radii from the central source, their 
result suggests a characteristic size that appears 
to decrease with increasing accretion rate or lumi- 
nosity. The clear evidence for a break in the soft 
band PDS of Akn 564 suggests a characteristic 
size associated with the soft X-ray excess emission 
which can be the size of the optically thick corona 
or the radius of the inner truncated disk in the 
model suggested by Done & Kubota (2006). It is 
also possible that the disk surface below an inner 
radius is ionized and the temperature of the in- 
ner disk surface is higher than that expected from 
the standard disk. The inner disk surface may 
then be responsible for the soft excess emission. 
The lack of a clear break in the hard band PDS 
of Akn 564 may suggest that power-law emission 
is probably not associated with a characteristic 
size and it arises from an extended corona sim- 
ilar to a patchy corona above an accretion disk. 
Future long observations of Akn 564 are required 
to clearly demonstrate the presence or absence of 
a break in the hard band PDS. 
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6. Conclusions 

We have presented temporal and spectral study 
of XMM-Newton observations of Akn 564 and 
Mrk 1044. The main results are as follows: 

1. The 0.2 - 2 kcV and 2-10 keV band X- 
ray emission from both the NLS1 galaxies 
are strongly correlated without any time de- 
lay. However, the variations in the 0.2 — 
0.5 keV band are found to lead those in the 
4-10 keV band by 1768 ± 122 s (90%) in 
Akn 564. We do not find any correlation be- 
tween the 0.2-0.3 kcV and 5-10 keV X-ray 
emission from Mrk 1044 at zero time delay. 
The two bands appear to anti-correlate at a 
positive lag of ~ 1000 s and likely correlated 
at a lag of ~ 3000 — 5000 s. Long observa- 
tions of Mrk 1044 are required to confirm the 
time delay. 

2. The full band power density spectrum of 
Akn 564 has a break at (1.2±0.3) x 10" 3 Hz, 
corresponding to a time scale of ~ 1000 s, 
similar to the delay between 0.2 — 0.5 and 
4 — 10 keV bands. It is not clear if there is 
any physical relation between the two time 
scales. There is a clear evidence of a break in 
the soft band PDS but the hard band does 
not show a clear evidence for the break. 

3. The soft (0.2 - 0.5 kcV) band PDS of 
Akn 564 is significantly steeper and weaker 
than the hard (3 - 10 keV) band PDS. This 
implies that the power-law component is 
more variable than the soft excess compo- 
nent on shorter time scales. The implication 
of this result on the Comptonization mod- 
els is that the variations in the power-law 
component are not only caused by the vari- 
ations in the seed photons but also by the 
variations in the hot electron population. 
Thus electron injection on short time scales 
(< 1000 s) are required. This could be evi- 
dence for magnetic flares thought to power 
the corona above the accretion disk. 

4. The soft excess emission from Mrk 1044 
is featureless in the EPIC-pn data and is 
well described by the smeared absorption or 
optically thick Comptonization. The RGS 
data show evidence for warm absorbers in 



Akn 564. The EPIC-pn spectrum of Akn 564 
is well described by a complex model consist- 
ing of optically thick thermal Comptoniza- 
tion in a low temperature (~ 0.15 kcV) 
plasma and a steep power-law, modified by 
two phase warm absorber medium and the 
Galactic absorption. The ionized reflection 
and smeared wind models fail to describe the 
data satisfactorily. 

5. The temporal and spectral characteristics of 
Akn 564 and Mrk 1044 are consistent with 
a two component corona - a compact low 
temperature, optically thick corona and an 
extended hot corona. The compact corona 
could be an inner optically thick region cou- 
pled to a truncated disk or the ionized sur- 
face of an inner untruncated disk. The soft 
excess emission can be produced in the inner, 
optically thick corona and the power-law 
component is produced in the hot corona. 
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